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I.  INTRODUCTION 

This  report  is  concerned  with  a star  identification  problem 
associated  with  star  mappers  on  a spacecraft  whose  orientation  is  not 
known.  Knowledge  of  spacecraft  attitude  can  be  lost  after  orbit  injec- 
tion or  after  unsupervised  or  uncontrolled  maneuvers.  In  any  case, 
the  basic  problem  is  that  little  or  no  a priori  information  exists  on 
the  inertial  pointing  direction  of  the  star  mappers  during  an  observa- 
tion period.  Under  these  "blind"  conditions,  it  is  not  clear  how  to 
associate  the  outputs  of  star  mappers  with  a catalog  of  known  stars. 

The  identification  problem  is  focused  on  star  mappers  which 
employ  a fixed  pattern  of  slits  and  allow  star  images  to  pass  to  a 
common  detector  (e.g.,  a photo  multiplier  tube).  In  this  case,  the  detec- 
tor emits  a series  of  pulses,  each  of  which  provides  a measure  of  the 
star's  visual  magnitude  and  slit-crossing  time  relative  to  a clock  (i.e.,  the 
transit  time).  However,  the  particular  slits  associated  with  the  pulses 
are  not  known.  Thus  the  outputs  of  these  star  mappers  have  an  inherent 
ambiguity  as  to  the  slits  associated  with  the  star  transits. 

The  report  presents  a star  identification  method  applicable  for 
star  mappers  with  slit  ambiguity  under  unknown  attitude  conditions. 

The  approach  assumes  that  the  spacecraft  can  be  controlled 
to  spin  about  a known  body  axis  prior  to  the  observation  period.  When 
this  condition  is  achieved,  star  mapper  outputs  are  obtained  and  pro- 
cessed to  provide  slit  and  star  identification.  Once  identified,  the 
transit  data  are  processed  for  attitude  determination. 

Star  identification  and  attitude  determination  were  undertaken 
by  Kenimer  and  Walsh^  for  a spacecraft  undergoing  torque-free 
spinning  motion.  Grosh^  generalized  their  results  by  allowing  for 
spacecraft  precession.  Both  works  employed  star  mappers  whose 
slits  allow  star  images  to  individual  detectors,  resulting  in  no  ambiguity 
as  to  the  particular  slit  generating  the  output  pulse.  The  proposed  method 
generalizes  their  approach  by  allowing  for  slit  ambiguity  and  by  accounting 
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for  the  possibility  of  incorrect  identification  due  to  measurement  and 
modeling  errors. 

The  development  begins  with  a description  of  the  star  mapper 
configuration  under  consideration.  Coordinate  systems  are  defined 
to  specify  the  orientation  of  the  spacecraft  and  star  mappers.  Relations 
describing  the  star  crossing  geometry  of  the  star  mappers  are  developed 
and  specialized  to  predict  parameters  associated  with  the  transit  data. 
The  statistics  of  these  parameters  due  to  measurement  and  modeling 
errors  are  developed.  These  relations  are  used  to  form  the  basis  for 
star  identification  and  gross  attitude  determination.  A brief  numerical 
study  is  presented  in  the  final  section. 


U.  SYSTEM  DESCRIPTION  AND  ASSUMPTIONS 


The  basic  elements  of  the  star  mapper  include  an  optical 
system,  reticle  and  a photoelectric  detector.  As  shown  in  Figure  1, 
the  optical  system  images  a field  of  stars  onto  a reticle  at  the  focal 
plane.  The  reticle  has  a fixed  pattern  of  slits  that  allows  star  images 
to  pass  to  the  detector.  As  star  images  are  carried  across  the  slit 
pattern,  the  detector  emits  a series  of  pulses  to  a processor  which 
computes  the  transit  time  and  visual  magnitude  corresponding  to 
each  pulse. 


The  particular  pattern  of  slits  under  consideration  is  shown 
in  Figure  2.  This  pattern  is  distinguished  by  the  presence  of 
two  parallel  slits  and  at  least  one  non-parallel  slit.  In  this  case,  the 
parallel  slits  are  designated  as  slits  1 and  2.  To  avoid  the  ambiguity 
between  the  slits  in  a region  about  the  optical  axis,  it  is  assumed 
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Figure  2.  Slit  Pattern  of  Sensor  j (Optical  Axis  View) 


It  is  assumed  that  there  are  J star  mappers  on  the  spacecraft. 

In  this  case,  a particular  star  mapper  is  referred  to  as  sensor  j. 
where  j = 1,  2,  • • •,  J.  For  simplicity,  each  star  mapper  is  presumed 
to  have  the  slit  pattern  of  Figure  2,  where  an  individual  slit  is  denoted 
as  slit  i (i  = 1,  2,  3,  4). 

Star  measurements  are  assumed  to  be  obtained  from  a spacecraft 
nominally  spinning  at  a constant  rate  u)q  about  a known  body  axis.  Each 
star  mapper  has  a fixed,  known  orientation  with  respect  to  the  space- 
craft. Deviations  from  the  nominal  spin  motion  and  star  mapper  geometry 

* The  development  can  be  generalized  to  incorporate  different  slit 
patterns  for  each  sensor. 


§ii  -8- 


are  assumed  to  be  small  over  a spin  period.  Denoted  as  6^,  these 
modeling  errors  are  assumed  to  be  independent,  normally  distributed 
random  variables  with  zero  means  and  covariance  P^. 


Sensor  j is  assumed  to  detect  stars  of  magnitude  or 

brighter  to  produce  a series  of  transit  time  measurements 

(k)  (k) 

t.  . and  visual  magnitude  measurements  M:.'  where  k denotes  the 

lJ  XJ  (k) 

index  of  the  observed  star.  The  errors  in  t. . ' are  assumed  to  be 

ij 

independent,  normally  distributed  random  variables  with  zero  means 

(k) 

and  covariance  P . The  errors  in  M:^  are  assumed  to  be  independent, 

uniformly  distributed  random  variables  with  zero  means  and  an  error 

(k) 

range  AM  dependent  on  M 


III.  COORDINATE  SYSTEMS 


For  the  assumed  spacecraft  motion,  it  is  convenient  to  specify 
the  orientation  of  the  spacecraft  and  star  mappers  in  terms  of  the 
following  coordinate  frames: 

• Earth-centered  inertial  reference  frames  (I  frame) 

• Spacecraft-centered  spin  reference  frame  (O  frame) 

• Spacecraft-centered  body  reference  frame  (B  frame) 

• Sensor-centered  boresight  reference  frame  (S^  frame) 

• Sensor-centered  trajectory  reference  frame  (N^  frame) 

• Sensor-centered  slit  reference  frame  (L^  frame) 

The  spacecraft  attitude  is  specified  by  the  I,  O and  B frames  whose 
basis  vectors  are  the  set  of  orthogonal  unit  vectors  (Xj,  y^,  z^), 

(5q.  yQ»  zQ)  and  (*B>  Yb*  resPectively*  As  indicated  in  Figure 

3,  the  unit  vector  z is  directed  toward  Aries,  Yj  is  along  the  North 
Celestial  Pole  and  x is  in  the  Celestial  Equatorial  plane.  The  unit 

A ^ A A t 

vector  Zq  is  along  the  nominal  spin  direction  and  xq  and  Yq  are  in  the 
spin  plane.  The  transformation  matrix  from  the  I frame  to  the  O frame 
is  denoted  as  co/i  and  can  be  expressed  as 

Cic»/I  = ^es^z  ^6s-*x  ^s-ly 


Figure  3.  Coordinate  System  Geometry  (I  and  O Frames) 


where  6 , £i_,  and  0„  denote  known  body  angles.  The  above  definition 

B 13  B 

allows  flexibility  in  specifying  the  spin  direction  relative  to  the  B frame. 

It  is  desirable  to  describe  the  orientation  of  sensor  j in  terms  of 
two  coordinate  systems:  the  S^.  frame  and  the  hh  frame.  The  basis 

vectors  of  the  S.  and  hT  frames  are  denoted  as  (xg  j,  ySj*  Zgj)  and 
(xN  .,  yN .,  z^),  respectively,  in  which  xg^  and  x^  are  along  the 
boresight  of  sensor  j and  yg^,  Zgj.  YNj,  are  in  the  focal  plane 

of  sensor  j.  The  frame  is  specified  relative  to  the  B frame  by  the 
transformation  matrix  Cg.  given  as 


CSj/B  = ["Pj]y  [Yj]z 


where  B and  y.  denote  boresight  angles  of  sensor  j.  On  the  other  hand, 

j j 

the  N frame  is  specified  relative  to  the  O frame  by  the  transformation 

j 

matrix  C^jj/O  &*ven  as 

CNj/0  = [_Ebj^y  [Abj]z  (?1 


where  Eb.  and  denote  boresight  elevation  and  azimuth  angles  of 

sensor  j.J  Since  = *Nj,  the  angles  Ebj  and  Abj  can  be  related  to 

0-,  V •,  6^,  and  0^  as 
J J 


■ 


■!. '' 


IV.  STAR  CROSSING  GEOMETRY 

As  the  spacecraft  spins,  the  trajectory  of  a star  may  cross 
a number  of  slits  on  the  star  mapper.  For  the  assumed  slit  pattern, 
this  number  ranges  from  one  to  four  depending  on  the  spin  direction 
and  the  sensor's  orientation  relative  to  the  B frame.  In  particular, 
there  are  two  basic  slit  crossing  geometries  which  are  of  practical 
inte  rest: 

Geometry  I:  3 Slit  Crossings/Star  Trajectory 

Geometry  II:  4 Slit  Crossings/Star  Trajectory 

Loosely  speaking,  one  of  the  two  non-parallel  slits  is  nearly  parallel 
to  the  y^j  direction  in  geometry  I;  none  of  the  slits  is  oriented  along 
y^.  in  geometry  U.  Since  a star  trajectory  in  the  y^  - plane  is 
approximately  parallel  to  y^,  the  majority  of  the  star  trajectories 
crosses  three  slits  in  geometry  I and  four  slits  in  geometry  II. 

For  a given  spin  direction,  those  sensors  resulting  in  geometry 
I or  II  are  the  sensors  under  consideration  in  the  remaining  develop- 
ment. For  these  sensors,  the  following  transit  time  difference  can 
be  predicted  from  the  star  crossing  geometry: 


A,  - ItW  t(k) 

AtDj  ' U3j  " l4j 


At 


Ej 


1 - t{$\.  geometry  I 

1 4^  " 4?  I • seometry  11 


>15) 


(16) 


* There  is  the  possibility  of  a two  slit  crossing  geometry  where  the 
two  parallel  slits  are  nearly  parallel  to  y^..  However,  this  geometry 
does  not  permit  any  means  to  discriminate  between  slit  crossings 
above  and  below  the  x^  - y^  plane. 

**  In  actuality,  the  nominal  star  trajectory  in  the  yNj  - z^  plane  is  a 
curve  with  mirror  symmetry  about  the  z^  axis. 
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where  At^  denotes  the  time  difference  between  the  parallel  slits 
and  At^  denotes  the  time  difference  between  the  extreme  parallel 
and  non-parallel  slits  of  a star  trajectory.  Expressions  used  to 
predict  Atjy.  and  At^  are  developed  in  this  section. 


(k)  (k)  /M 

The  terms  E:.'  and  A.,  denote  elevation  and  azimuth  angles  of  uv  ' 
ij  ij  6 

defined  relative  to  the  O frame  (see  Figure  7).  The  elements 
of  are  specified  as 

qOij  ° C0/B  CS;/B  CLij/Sj  qLij 


qur  ° 


Substituting  (19)  - (22)  into  (17)  yields 

(k)  „(k)  .(k)  (k)  _(k)  . A(k)  (k)  . _(k) 

q . . cos  E. . cos  A. . + q . . cos  E. . sin  A. . + q . . sin  E. . 

xij  ij  IJ  >lj  ij  IJ  MZlj  ij 


Through  the  relations  of  (5),  (6)  and  (14),  it  is  seen  that  (25) 
(k)  (k) 

involves  E; . , A..  , 6D,  9Dt  3-»  Y-»  6..  and  o...  Among  these 

ij  ^ B B B J J U !J  /kv  /k) 

parameters,  on,  Qn,  0_,  3 . and  y . are  known  and  E. . , A..,  6.. 
r B B’  B’  j j ij  ij  ij 

and  a.,  are  unknown, 
ij 


Figure  7.  Star  Crossing  Geometry 
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B.  SLIT  RELATIONS 

In  general,  the  elevation  and  azimuth  terms  in  (25)  can  be 
expressed  as 

E*k)  = E + T]!k) 

bj 


A<k) 


A . + 
bj  ij 


where  T)!k^  denotes  the  elevation  of  p/k^  relative  to  the  x,. . - y,.. 

# ^ INJ  INJ. 

plane  and  £>..  denotes  the  azimuth  angle  defined  relative  to  A^;  m 


ij 


the  Xq  - y q plane.  Also  it  is  noted  that 


T ^ r T 

CNj/0  °Nj/Sj  ^Lij/Sj  qLij 


Substituting  (26)  - (28)  into  (25)  yields 


dxi  COS  <£bj+T1Sf)  COS  i^yi  C°'<Ebj+10  Sin  Ai] 


+ d . sin  (E,  . +T ]!k))  = 0 
zi  bj  ‘ij 


(k) 


whe  re 
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rr;  *r'  r*  nr  r-«r . W 


(26) 

(27) 


(28) 


(29) 


d . 

XI 

= cos  E,  . sin  6. . - 
bj  ij 

sin  E,  . 
bj 

COS  t . . 

IJ 

cos  6 . . 
ij 

(30) 

• 

d . 
Y1 

= - sin  . . cos  6. . 

VJ  ij 

(31) 

* . 

dzi 

= sin  E,  . sin  6. . + 
bj  ij 

cos  E,  . 
bj 

cos  t ■ • 
IJ 

cos  6. . 
U 

(32) 

y 

= Of.  . - p . 

ij  J 

(33) 

Furthermore,  7 y * can  be  expressed  as 


= arc  sin  [sin  t,:  sin  cp!.  - cos  i|(.  . sin  6..  cos  cp'kH 

*J  ij  ij  Yij  J 


where  the  terms  cpjk  denotes  the  angle  between  and  xT  . 

ij  Y Lij 

Suppose  star  k crossed  slit  i at  time  tf1^  and  slit  i'  at  time 

t|!^.  The  predicted  time  difference  between  t!!^  and  tfk^  can  be 

J 1 J U 

expressed  as 

(k,  ^ - A<k> 

At(k)  =-^J LL 

l'/l  (U 


Given  slit  i andcp:.  , 

XJ  (k) 

one  can  first  determine  J\)  ' using  (34)  and  then  numerically  solve 
(k)  « 3 

(29)  for  A\  . Similarly  with  i = i',  one  can  numerically  solve  (29)  for 

A.t.  under  the  constraint  T|j,)  = T)|.^  and  subsequently  determine  At^) 

• „rl  1 J i'/i 

using  (35). 


(k)  (k) 

Because  cp^ j or  T]^.  is  not  known,  the  predicted  time  diffe 


rences 


f°r  AtDj  and  AtEj  are  aPPr°ximated.  In  particular,  the  time  At^  between  the 
parallel  slits  is  predicted  from  (35)  under  the  condition 


"Si?  - = ° 


(i1,  i = 3,  4) 


The  resulting  value  of  AtDj  represents  an  average  over  the  field  of 
view  limits. 

The  time  At£^  between  the  extreme  slits  is  predicted  from 
(35)  under  the  condition 

(k) 

*ij  =-FU  s®n  (£bjJ 
T]lk)  = n<k> 

1 j ij 


1 


where 


i'.  i = 


2,  4 

1.  4 


geometry  I 
geometry  II 


(39) 


and  sgn  (x)  denotes  the  sign  of  x.  Thus  At^  and  At£^  can  be  predicted 


directly  from  the  star  crossing  geometry. 


ARC  LENGTH  RELATIONS 


q<k)  = q<k> 
qi'  qOi' 


-cos  9^k)  sin 

i'i  i'i 


,i»  »!|° 

i'i 


cos  0<k) 

l'l 


A similar  group  of  expressions  can  be  obtained  at  times  t!f^  and 

I Qj\  1 J 1 0 \ 

t.,:  which  corresponds  to  star  l.  To  be  specific,  the  terms  £:.  , 

•Mi)  lJ 

, A\  and  A^  are  described  by  (40)  - (48)  with  i substituted 

for  k. 

The  arc  length  between  stars  k and  i is  defined  as 


= arc  cos 


[alk»  • sw] 


The  O frame  components  of  at  time  t^  and  at  time 

ij  ij 

have  the  form  of  (19).  Substituting  these  components  into  (49)  yields 


d,  . = arc  cos  sin  E^  sin  e!^ 
k*  n n 


+ cos  e!^  cos  e!^  cos 


[A!Jk,-A!j,4”oC-tSj,il!  (50> 


— (k)  (k)  (/) 

Thus  d,  . is  expressed  in  terms  of  the  transit  times  t. . , t., . , t. . and 

(i)  * lJ  1 J 1J 

t:,..  In  contrast,  the  arc  length  of  stars  o and  P from  a star  catalog  is 

expressed  as 


* Equation  (50)  is  readily  extended  to  the  situation  where  the  transits 
result  from  two  different  sensors. 


d „ = arc  cos  }sin  6 sin  6Q 
ag  ) a;  P 


D.  ERROR  RELATIONS 

The  previous  sections  indicate  that  Atj-j.  and  d^  can  be  expressed 
as 


1 


MDj  = AtDj  (0'  T1)  (52) 

dkX  = <@*  T>  <53) 

where 


(54) 

(55) 

(56) 

If  ©^  and  Tlj^j  denotes  the  nominal  values  of  © and  T|,  respectively, 
deviations  about  the  nominal  can  be  expressed  as 


0 [6ij'  6irj'  “ij*  “i'j*  Pj * Yj’  V nP'  V 


71  = t^if'  TliJj  3t 


- [t(k),  e>,  t<x>,  tift 

lj  l'j  1J  l'j 


6©  = 0 - 0N 


(57) 


ST]  = T1  - 11N 

where 

nN  = [o,  o]T 


(59) 

(59) 


The  terms  60  and  6T)  represent  system  model  errors.  Measured  transit 
times  are  represented  by  t and  their  errors  are  denoted  as  6t. 

To  approximate  the  effects  of  model  and  measurement  errors 
on  and  dk£*  (52)  and  (53)  are  linearized  about  ©j^,  and  t to 

yield 


“Dj  ' J«9  60  + JtT|  6” 


d.  , = J 6©  + J,  6t 

kz  d®  dT 


(60) 

(61) 
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where 


° &tDj 


t©  <)  © 
d Ot 


tt]  a t\ 

d d 


2i 


k4 


d©  <3  © 

3 d 


©N*  \ 


®N* 


®N*  T 


W 


dT  9 T | @N>  T 


(62) 


(63) 


(64) 


(65) 


In  this  development,  60,  6 Tj  and  6t  are  assumed  to  be  normally  distri- 
buted random  variables  with  zero  means  and  the  following  second-order 
statistics: 


T 

E [6©  6©  ] 

® 

Oh 

II 

(66) 

E [611  5T]T] 

~ PT1 

(67) 

E [6t  5t 

= P 

T 

(68) 

E [ 6T)  6@T ] 

= 0 

(69) 

E [ 6T|  6tT] 

= 0 

(70) 

E [6t  6©T] 

= 0 

(71) 

where  P„P  and  P denote  covariance  matrices  of  6®,  6t  and  6T1, 

fcT  T T| 

respectively.  The  error  covariances  P^  and  P^  are  specified  by  the 
problem  at  hand.  The  error  covariance  P^  can  be  modeled  as 


A1 


2x2 


(72) 
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when 


°n  = <FU>-  Vj  (FU>’ 


The  terms  T|„  and  are  computed  from  (34)  with  <p„  = cp^  = F^; 


I is  an  nxn  identity  matrix, 
nxn 


Based  on  the  first-order  approximations  of  (60)  - (61)  and 
the  statistics  of  (66)  - (71),  A and  d,  , become  normally  distributed 

with  zero  means  and  variances  ct^  andCT^^,  respectively,  where 


ctD  = Jt©  P©  J©  + JtT|  PT)  JtTl 


2 T T 

of#  = JT0+  J,  P JT 

k4  d®  © d®  dT  t dr 


rrarr\  « 


The  approach  to  star  identification  begins  with  a set  of  star 
measurements  obtained  from  a spinning  spacecraft.  In  this  case,  the 
relative  times  between  successive  slit  crossings  are  predicted 
from  the  star  crossing  geometry.  These  predictions  are  compared 
to  transit  time  data  to  provide  a basis  for  slit  identification.  Those 
data  matching  the  time  predictions  and  having  similar  visual  magni- 
tudes are  designated  as  "multiple"  transit  data.  The  remaining 
data  are  designated  as  "single"  transit  data. 

Once  the  slits  are  identified,  the  arc  lengths  between  stars 
are  computed  from  the  multiple  transit  data  and  compared  to  those 
from  a star  catalog.  Among  the  star-pair  combinations  matching 
the  data,  a selection  is  made  on  those  stars  that  form  a consistent 
star-pair  chain. 

A flow  chart  of  the  star  identification  method  is  shown  in 
Figure  8 . It  is  noted  that  the  method  consists  of  three  basic  opera- 
tions: (1)  an  initialization  operation  which  inputs  measurement  and 
model  data  and  computes  system  parameters,  (2)  a screen  operation 
which  identifies  and  isolates  multiple  transit  data  and  (3)  a star 
selection  operation  which  identifies  stars  corresponding  to  the  multiple 
transit  data. 


A.  INITIALIZATION 

The  first  task  is  to  input  system  data,  a star  catalog  and  star 

mapper  data.  System  data  include  star  mapper  geometry  parameters 

(Vj,  a.j,  6^,  Fy,  F L),  spin  motion  parameters  (ujq,  6^,  Q^,  0g), 

and  error  parameters  (Pg,  PT»  AM).  The  star  catalog  consists  of 

line  of  sight  parameters  (6  , Q ) and  the  cataloged  visual  magnitude 
(a)  a a 

(M^  ')  associated  with  star  a.  Once  obtained,  the  star  catalog  is 
reduced  by  accounting  for  the  maximum  sensitivity  of  the  star  mappers, 
(i.e.  , M*k)  £ Mcm)»  any  a priori  knowledge  of  the  inertial  spin  direc- 
tion, and  the  effects  of  star  obscuration  (e.g.,  by  the  Earth).  The 
reduction  results  in  a limited  star  catalog  which  is  used  for  the  star 
selection  task. 

The  star  mapper  data  consists  of  a series  of  transit  times 
(k)  (k) 

and  visual  magnitudes  (t^  , M:/).  For  this  discussion,  it  is  assumed 
that  the  raw  data  have  been  processed  over  a number  of  spin  periods 
(scans)  to  reduce  the  amount  of  spurious  transits.  To  this  end,  multi- 
scan techniques  can  be  used  if  the  spacecraft  does  not  appreciably 

(3) 

precess  during  successive  scans.  Once  completed,  the  resulting 
data  are  reduced  to  a single  scan  for  star  identification. 


V 1 
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B. 


SCREEN  OPERATION 


The  screen  operation  consists  of  a series  of  tests  to  isolate 
and  identify  multiple  transit  data.  The  first  test  is  aimed  at  isolating 
those  pairs  of  transits  associated  with  a star  crossing  the  parallel 
slits  (i.e.  , slits  3 and  4).  If  two  transits  are  from  the  same  sensor 
and  within  a f ield-of-view  time  window,  then  they  qualify  as  a candi- 
date pair.  Given  a candidate  pair  with  transit  times  t„.  and  t . and 

mj 

magnitudes  M.  . and  M from  sensor  j,  a time  residual  r^.  and  a 

mj  Dj 

magnitude  residual  r^  are  computed  as  follows: 


rn  = t . - t . . - At_  . 
Dj  mj  Aj  Dj 


(76) 


r.  . . = M.  . - M. 
Mj  mj  Aj 


(77) 


where  At^j.  is  specified  by  (35)  and  (36).  Also  a time  score  and 
a magnitude  score  are  computed  as 


_ 2 .,-1 
CDj  " rDj  Vt 


(78) 


CMj  I rMj  ’ 


(79) 


where 


it  _ o 2 2 

Vt  - 2 CTt  + CTtD 


(80) 


The  terma  is  a specified  variance  of  the  transit  time  measurement  error  and 
2 1 

atD  *s  comPuted  from  (74). 

The  decision  rule  is  to  accept  the  times  (t.  t .)  as  "double 

mj 

slit"  times  if  the  following  inequalities  are  true: 


CDj<  CDB 


(81) 
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n ztj. 


(82) 


where  is  a threshold  level  for  the  double  slit  test.  Because  of 

the  assumptions  on  the  error  statistics,  the  residuals  r^  is  normally 
distributed.  Thus  the  score  has  a chi-squared  distribution  with 
one  degree  of  freedom.  Moreover,  given  a specified  probability  of  a 

correct  decision,  the  threshold  can  be  obtained  from  standard 

(4)  DB 

tables.'  ' 


The  next  task  is  to  isolate  the  remaining  transits  associated 
with  the  double  slit  transits.  To  this  end.  At,-.,  is  computed  from  (35)  and 
(37)  - (39)  to  provide  a time  window.  Transit  times  are  checked  on  both 
sides  of  the  double  slit  times  to  determine  if  they  are  within  the 
window.  Those  candidates  within  the  window  are  tested  for  similar 
visual  magnitudes  as  in  the  double  slit  test.  Furthermore,  the  number 
of  candidates  within  the  window  are  checked  for  consistency  with  the 
known  star  crossing  geometry.  The  allowable  number  of  candidates/ 
windows  is  one  in  geometry  I and  two  in  geometry  II.  If  there  is  an 
allowable  number  on  both  sides  of  the  double  transits,  both  candidates 
are  rejected  to  avoid  ambiguity. 


Those  candidates  with  the  allowable  number  are  combined 

with  the  double  transits  to  form  the  k—  multiple  transit  group.  Slit 

(k)  (k) 

designations  are  assigned  and  the  parameters  A.,  and  E..  are 
computed  via  (43)  - (48).  The  above  process  continues  until  all 
of  the  transits  have  been  considered. 
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c. 


STAR  SELECTION  OPERATION 


The  star  selection  operation  consists  o£  an  arc  length  test, 
a star-pair  consistency  test  and  a magnitude  test. 

The  arc  length  test  begins  by  computing  the  arc  length  between 
two  stars  associated  with  the  first  two  multiple  transit  groups.  In 
general,  for  groups  k and  (k  / l) , the  arc  length  d^  is  computed 
via  (50)*  Then  a residual  r^  is  computed  from  d^  and  the  arc  length 
d - between  stars  a and  6 of  the  limited  star  catalog  as 

aP 


(83) 


where  d Q is  specified  in  (51).  Further- 

O'  P Q 

more,  a chi- squared  score  based  on  the  residual  is  computed  as 


/ OfP. 
(rki> 


Pv2 


(84) 


where  is  specified  in  (75). 

The  decision  rule  is  to  accept  the  stars  a and  P as  candidates 
if 

< CAB  <85> 


where  is  a threshold  selected  to  meet  a specified  probability  of 

correct  decision;  otherwise,  the  star-pair  is  rejected.  Every  star- 
pair  combination  in  the  limited  catalog  is  tested  in  this  manner. 

Once  the  star  candidates  for  groups  (k,  l)  are  obtained,  the 
same  procedure  is  used  to  obtain  candidates  for  groups  (£,  m)  and 
(k,  m).  Given  star  candidates  for  three  groups  of  data,  a search 
is  conducted  to  identify  those  candidates  which  are  logically  consistent 
with  the  data.  For  example,  stars  a,  P and  y comprise  a consistent 
triplet  if  each  is  a candidate  in  only  two  of  the  groups  (k,  l),  (l,  m) 
and  (k,  m). 
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Finally,  each  of  the  candidates  comprising  a consistent  triplet 

is  subjected  to  a magnitude  test.  In  this  case,  a residual  is  computed 

from  the  cataloged  visual  magnitude  of  star  a and  the  average 

tli  ^ 

magnitude  associated  with  the  k — multiple  transit  group,  or 

= M<“>  - Mlk)  (86) 


where 


n, 


M 


(k)  - f M<k)/n, 
i = 1 1J  k 


(87) 


and  n^  denotes  the  total  number  of  transits  for  the  k—  group.  The 
rule  is  to  accept  the  three  candidates  as  stars  associated  with  the 
multiple  transit  group  k if  all  three  stars  satisfy 


l*Mkl  < AM 


(88) 


In  a similar  manner,  star  triplets  are  identified  for  the  next 
three  transit  groups  (i.e.,  the  last  two  of  the  previous  group  and  one 
additional  set).  The  process  continues  until  all  of  the  transit  groups 
are  considered.  Each  candidate  star  is  thus  required  to  form  a con- 
sistent triplet  with  three  other  pairs  of  candidates- -the  preceding  pair, 
the  subsequent  pair,  and  the  pair  just  before  and  after  the  candidate. 
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VI.  ATTITUDE  DETERMINATION 


The  problem  is  to  determine  an  estimate  of  the  spacecraft's 
inertial  attitude  based  on  star  mapper  measurements.  It  is  noted 
that  the  transformation  matrix  between  the  B and  I frames  can  be 
expressed  as 


CB/I  “ C0/B  *0B'  UB’  6B*  C0/I  (0s’  Qs’  6 s* 


(89) 


Assuming  that  0B>  and  of  are  known,  the  problem  is 


reduced  to  determining  6g,  Qg  and  9? 


A gross  estimate  can  be  obtained  by  utilizing  line  of  sight 
information  of  stars  identified  with  the  multiple  transit  data.  Speci- 


-(a) 


ficaily,  the  inertial  components  of  p can  be  expressed  as 


» 


'O/I 


(a) 


(6s*  Qs’  es) 


(90) 


where 


(a)  _ 


cos  6 sin  Q 
a a 


sin  n 


a 


cos  6 cos 

a a 


(91) 


and  is  specified  in  (21)  in  terms  of  a|®'  and  which  are  known  for 

each  multiple  transit  group.  Thus  (90)  provides  two  independent  relations 


for  each  identified  star. 


The  undetermined  parameters  6g,  Ug  and  0g  can  be  estimated 


based  on  the  following  model: 
6 (t)  = 6 


so 


(92) 


a (t)  = a 

s so 


0 (t)  = a)  (t  — t ) +9 
s o'  o so 


(93) 

(94) 
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where  5 , Q and  9 are  constants  corresponding  to  time  t . 

so  so  so  r a Q 

Thus  given  N identified  stars,  the  estimation  problem  can  be 

formulated  to  select  parameters  6 , 0 and  6 such  that  thev 

so  so  so  7 

minimize  the  estimation  error 


e = E (pm  " ?rr) 
m = 1 

subject  to  the  constraints  of  (90)  - (94)  and 


The  term  p denotes  the  estimate  of  p where  m = 1,  2,  3.  Given 

at  least  two  identified  stars,  the  estimates  can  be  computed  via  an 

iterative  least- squares  procedure.  Once  pm  are  known,  estimates 

of  6 , w and  0 are  computed  as 
s s s 


6 = 6 
s so 


A A 

q = n 

s so 


9 =o)  (t-t)  + e 

« o o so 


Finally,  the  estimate  of  the  spacecraft's  inertial  attitude  is 
represented  by  the  matrix 


CB/I  = C0/B  C0/I  (9s'  ns'  6s> 


A reasonably  close  initial  estimate  of  the  unknown  parameters 

can  be  obtained  from  the  data  of  two  identified  stars.  Given  M.n^» 

la)  (01  ^ 

Pj  ' and  p^  ' of  stars  a and  3,  the  unknown  parameters  p^  and  p^  can  be 

computed  as  follows: 

pt  = arc  sin  (s^)  (101) 

p2  = arc  tan  (sIx/sl2)  (102> 

where  s^,  s^,  Sjz  are  the  I-frame  components  of  the  spin  direction, 


SI  = (sIx’  Sly'  Slzl 


(103) 


But  Sj  can  be  expressed  as 


5i _ n '"Si  sSl 


(104) 


where 


sSI  = (0  , 0 , 1)T 

(105) 

a _ r (o)  (3)  n(a3)n 

ASI  ~ L^SI  ^SI  SI  J 3x3 

(106) 

A = \Ja)  u(P)  n(o0)l 
Aj  Lm-j  Pj  nj  J3x3 

(107) 

Pc**  = O t ) + 9 ]T 

i>I  O 1J  o soJz 

(a) 

M-  o 

(108) 

pi^  = [u>  (t!P)  - t ) + 9 JT 

^SI  O lj  o so J z 

JP) 

Po 

(109) 

The  terms  nj.j^  and  nj°^  are  inertial  components  of  the 

A(arf3)  , . , . ..  , . ~(a)  , ->(3) 

vector  n which  is  orthogonal  to  p and  p , or 


(110) 


'(<*0)  - u(g)x  u(g) 

n l^x^l 


Thus  for  a given  0 the  elements  of  u,g^  and  and  hence 

and  can  be  computed.  In  turn,  an  estimate  of  p^  can  be  determined 

as 


p^  = arc  tan 


(a)  (ot)  (o) 

sin  p1  sin  p2  p}x  + cos  Pj  p,,  + sin  Pj  cos  p£  p,^' 


cos  p2  ^Ix  ‘ sm  P2  »lz 


Ja)  - oo  (t<a)  - t ) 

1J  O lj  o 


(111) 


Once  an  attitude  estimate  is  obtained,  a detailed  star  identifi- 
cation operation  can  begin.  A typical  procedure  is  to  predict  transit 

times  based  on  the  gross  attitude  estimate  and  to  match  the  results 

(3) 

with  the  measured  single  transit  data. 
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VII.  NUMERICAL  STUDY 


A brief  numerical  study  was  conducted  to  test  the  feasibility 
of  the  proposed  star  identification  method.  To  this  end,  a computer 
program  with  the  basic  operations  was  implemented  and  tested  with 
simulated  star  data.  The  tests  consisted  of  generating  star  mapper 
data  under  known  conditions  and  then  processing  the  data  through  the 
program  to  identify  stars  and  estimate  spacecraft  attitude.  A 
measure  of  performance  was  obtained  by  comparing  the  program's 
results  to  the  true  star  observations  and  attitude  motion. 
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A.  RUN  CONDITIONS 

Star  data  were  generated  for  two  basic  types  of  spacecraft 
motion:  (1)  pure  spin  motion  (designated  as  motion  A)  and  (2)  nuta- 
tion motion  (designated  as  motion  B).  In  motion  A,  the  spin  rate 

was  set  to  id  + Au>  , where  Au)  denotes  a spin  rate  error  assumed 
o o o 


— = 2n  x 10  ^ 


Thus  the  actual  spin  angle  was  simulated  as 


9 =9  + (u)  + Au>  ) (t  - t ) 

s so  o o o 


In  motion  B,  spacecraft  was  assumed  to  cone  and  precess 

away  from  the  nominal  spin  direction.  In  particular,  the  half-cone 

angle  was  set  to  |6,.|  and  the  cone  frequency  was  2 cd  ; the  precession 
“4  ^ 

rate  was  2x10  cd  . 

o 

Three  star  mappers  were  assumed  to  be  oriented  on  the  space- 
craft such  that  each  sensor  viewed  a different  portion  of  the  celestial 
sphere  during  a scan.  The  spin  direction  in  the  B-frame  was  selected 
to  result  in  slit  crossing  geometry  I for  the  three  sensors.  The  inertial 
spin  direction  was  assumed  to  be  directed  near  the  sun  line.  Thus, 
depending  on  the  day  of  the  year,  the  three  sensors  observe  a pre- 
dictable number  of  stars.  Assuming  that  each  sensor  is  sensitive  to 
4th  magnitude  stars  and  brighter  (i.e.,  Mcm  = 4),  Figure  9 indicates 
the  daily  number  of  stars  within  the  field  of  view  of  the  three  sensors. 


To  provide  a range  of  observed  stars,  star  data  were  generated 
for  January  12  (dense  star  day)  and  for  March  7 (sparse  star  day).  In 
both  cases,  data  were  based  on  the  following  models  of  the  measure- 


ment errors: 


9JE)s  9jqeT?*AV  J°  Jaquinfj  jbjox 
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Figure  9.  Star  Availability  versus  Day  of  Year 


(114) 


&M  = (.01  Mfk)  + .04)  M 

ij  i,ni 

at  = 6*10-3Ts 


where 


M 


cm 


= 4 


T = 2 16-,1/u) 
s 3j 1 o 


(115) 


(116) 

(117) 


Parameters  used  by  the  star  identification  program  were 
mismodeled  from  the  simulated  data.  In  particular,  the  program 
modeled  the  spin  rate  as  and  the  error  covariances  PQ  and  Pt 
as 


(118) 


P 

T 


2 

a 


T 


I4x4 


(11 9> 


where 


= 0 


a 


2 

T 


+ 


«ir)2 


(120) 

(121) 


Given  the  star  data,  the  program  operated  with  and  without  the  visual 
magnitude  tests.  For  the  transit  time  tests,  the  decision  thresholds 
were  selected  to  achieve  a .99  probability  of  correct  decision  (i.e., 

cdb  = cab  = 6-635)- 


B. 


RESULTS 


The  results  of  the  star  identification  runs  are  summarized 
in  Table  I.  The  table  indicates  that  no  slits  or  stars  were  incorrectly 
identified  for  the  range  of  parameters  considered.  Some  stars  were 
not  identified  because  of  single  transits  and  densely- spaced  multiple 
transits  with  similar  visual  magnitudes.  Without  the  visual  magnitude 
tests,  the  program's  capability  to  discriminate  between  densely- spaced 
data  was  further  impaired.  Also  nutation  tended 
to  reduce  the  number  of  identified  stars. 

The  rejection  of  densely- spaced  transits  with  similar  magni- 
tudes occurred  during  the  screen  operation  to  avoid  multi-star  crossing 
ambiguities.  This  result  implies  that  there  is  a star  density  limit  above 
which  no  stars  will  be  identified  by  the  method.  This  limit  can  be  approxi- 
mated by  requiring  that  the  star  density  dg  be  no  more  than  1 star  per  area 
covered  by  the  field  of  view,  or 


d £ d 
s sm 


(122) 


(123) 

(124) 


and  kp  is  a geometry  shape  parameter. 

Based  on  stars  from  the  SAO  catalog,  the  star  density  in  stars/ 

square  degree  brighter  than  M are  plotted  against  the  visual  magni- 

c (5) 

tude  Mc  are  plotted  in  Figure  10.  ' Curves  are  given  for  the  average  density  as 

well  as  the  densities  at  the  extreme  galactic  latitude  bands.  Super- 
imposed are  levels  of  dgm  based  on  a square  shape  area  (kp  = 2/tt) 
for  various  field  of  view  angles  Fy.  It  is  noted  that  as  the  field  of 
view  increases,  the  screen  operation  without  visual  magnitude  discri- 
mination tends  to  reject  dim  star  data.  Thus  for  data  processing 
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Figure  10.  Star  Density  versus  Visual  Magnitude 


can  be  reduced  to 


purposes,  the  visual  magnitude  parameter 
be  consistent  with  the  level  of  Figure  10.  This  reduction  will  limit 
the  number  of  stars  considered  from  the  star  catalog  and  improve 
the  computational  efficiency  during  the  star  selection  operation. 

Figures  11  - 19  indicate  the  transit  data  of  the  January  12 
run  at  various  stages  of  the  star  identification  program.  The  run 
assumed  motion  A with  visual  magnitude  tests.  The  first  stage  occurs 
at  the  input  to  the  program;  the  second  stage  occurs  at  the  output 
of  the  screen  operation;  and  the  third  stage  occurs  at  the  output 
of  the  star  selection  operation.  It  is  noted  that  the  screen  operation 
tends  to  reject  densely- spaced  dim  star  data.  The  reason  is  due  to 
the  stringent  tests  designed  to  provide  a high  probability  of  true 
identification  at  the  expense  of  missing  valid  stars.  In  effect,  the 
star  identification  program  tended  to  miss  a valid  star  rather  than 
accept  a false  star. 

The  rss  attitude  estimation  error  resulting  from  a least 
squares  solution  of  the  January  12  run  is  shown  in  Figure  20.  The 
error  levels  reflect  the  attitude  deviations  from  the  nominal  motion 
that  were  not  modeled  by  the  estimation  program.  Improvements  can 
be  expected  by  formulating  the  estimation  problem  with  a model  that 
reflect  these  attitude  deviations. 
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Figure  13.  Transit  Data  Passing  the  Star  Selection  Operation  (Sensor  1) 
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Figure  14.  Input  Transit  Data  (Sensor  2) 
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Figure  16.  Transit  Data  Passing  the  Star  Selection  Operation  (Sensor  2) 
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VIII.  SUMMARY  AND  CONCLUSION 


The  paper  has  presented  a method  to  identify  stars  and  slits 
associated  with  the  outputs  of  star  mappers  whose  inertial  orientation 
is  not  known.  The  approach  assumed  that  star  data  are  obtained  as 
the  spacecraft  performs  a spinning  maneuver  about  a known  body 
direction  but  an  unknown  inertial  direction.  Given  a known  orienta- 
tion of  the  star  mappers  relative  to  the  spacecraft's  body  frame, 
expressions  were  developed  to  predict  the  relative  transit  times 
between  successive  slit  crossings  and  to  describe  the  arc  length 
between  stars.  First-order  effects  of  modeling  and  measurement 
errors  on  the  resulting  expressions  were  obtained  to  provide  a statis- 
tical basis  for  slit  and  star  identification.  Also  a least-squares  pro- 
cedure was  outlined  to  estimate  the  attitude  of  the  spacecraft  based 
on  the  identified  star  data. 

A preliminary  numerical  study  indicated  that  the  star  identi- 
fication method  is  feasible  for  a limited  range  of  star  density.  Because 
high  star  densities  increase  the  possibility  of  multi-star  crossing 
ambiguities,  the  method  tended  to  reject  dim  stars  from  considera- 
tion. However,  those  stars  that  were  processed  were  correctly 
identified  for  the  range  of  parameter  considered. 


REFERENCES 


1.  R.  L.  Kenimer  and  T.  M.  Walsh,  "A  Star  Field  Mapping  System  for 
Determining  the  Attitude  of  a Spinning  Probe,"  Paper  presented  at 
the  Aerospace  Electro -Technology  Symposium  of  the  International 
Conference  and  Exhibit  on  Aerospace  Electro -Technology,  April  19- 
25,  1964. 

2.  C.  B.  Grosch,  "Orientation  of  a Rigid  Torque-Free  Body  by  Use  of 
Star  Transits,"  A LA  A Journal  of  Spacecraft  and  Rockets,  4^(5),  562- 
566  (May  1967). 

3.  C.  B.  Grosch,  A.  E.  La  Bonte,  and  B.  D.  Vannelli,  "The  SCNS  Atti- 
tude Determination  Experiment  on  ATS-III,"  Proceedings  of  the  Sym- 
posium on  Spacecraft  Attitude  Determination,  September  30  - 
October  1-2,  1969,  207-221  (31  October  1969). 

4.  Standard  Mathematical  Table,  14th  Edition,  ed.  S.  M.  Selby,  The 
Chemical  Rubber  Company,  Cleveland,  Ohio,  270-271  (1965). 

5.  J.  D.  Gilchrist,  Final  Report  on  Project  STRETCH  (Space  Technology 
Requirements  and  Engineering  Tests  of  Components  Hardware), 

Vol.  Ill,  TR-0066(57o3 )- 1 , The  Aerospace  Corporation,  El  Segundo , 
California  (July  31,  1970). 


